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Based on the exact relationship to Random Matrix Theory, we derive the probability distribution 
of the fc-th smallest Dirac operator eigenvalue in the microscopic finite- volume scaling regime of QCD 
and related gauge theories. 



PACS number(s): 12.38.Aw, 12.38.Lg 

While index theorems have for long been known to re- 
late the number of exact zero modes of the Dirac operator 
to gauge field topology, it is only during the 1990's that 
it has been realized how much further one can push such 
exact predictions. When the gauge theory in question 
is in a phase with spontaneously broken chiral symme- 
try, Goldstone's Theorem allows for much more detailed 
statements, beyond the zero modes. The bridge between 
the Dirac operator spectrum and the Goldstone degrees 
of freedom is the effective partition function of the asso- 
ciated chiral Lagrangian [1,2]. The infinite- volume chiral 
condensate will be denoted by E. By considering the the- 
ory at finite four- volume V , one can choose a microscopic 
scaling regime in which this volume is sent to infinity at 
a rate correlated with the chiral limit m — > such that 
fx = T,Vm is kept fixed. Then exact statements can be 
made for an infinite sequence of Dirac operator eigenval- 
ues that accumulate towards the origin. This realization 
was first made on the basis of a universal Random Ma- 
trix Theory reformulation of the problem [3-6] , but it has 
since then also been established directly at the level of 
the effective Lagrangian [7,8]. 

Conventionally, the exact analytical statements that 
can be made about this infrared region of the Dirac oper- 
ator spectrum in finite-volume gauge theories have been 
phrased in terms of the microscopic spectral n-point func- 
tions, and in particular the microscopic spectral density 
Ps{C',{p}) itself (here ( = EVA is the microscopically 
rescaled Dirac operator eigenvalue A) [3,4,9-18]. Also 
the exact distribution of just the lowest Dirac operator 
eigenvalue has been derived on the basis of the relation- 
ship to Random Matrix Theory [10,19-23,15,16]. In fact 
also the distribution of the second-smallest eigenvalue in 
what corresponds to the massless limit of SU(N C > 3) 
gauge theories with fermions in the fundamental repre- 
sentation can be extracted from the work of Forrester 
and Hughes [21]. In this paper, we push these compu- 
tations one step further and provide analytical expres- 
sions for both the joint probability distribution of the 
first fc eigenvalues and the distribution of the fc-th eigen- 



value. We do so for the chiral unitary ((3 — 2), orthogonal 
(}3 = 1, corresponding to SU{2) gauge theories with fun- 
damental fermions), and symplectic ensembles (/3 = 4, 
corresponding to SU(N C > 2) gauge theories with ad- 
joint fermions). We will treat the most general case of 
massive fermions, and (with one exception noted below) 
a sector of arbitrary topological charge v. In the micro- 
scopic scaling limit this can be done to arbitrarily high 
order fc as long as the volume V is taken large enough. 
This gives us an infinite sequence of distributions of Dirac 
operator eigenvalues that all, on the macroscopic scale, 
build up the spectral density at the origin p(0). We stress 
that there is much more information in these individual 
eigenvalue distributions than in the summed-over micro- 
scopic spectral density Ps(C; {/■*}) itself. In particular, 
from the point of view of lattice gauge theory simulations 
it is far more convenient to be able to perform compar- 
isons with individual eigenvalue distributions than with 
just the average eigenvalue density. 

There are by now detailed analytical predictions for 
the microscopic spectral densities of the Dirac operator 
for all three universality classes. These microscopic spec- 
tral densities show a typical oscillatory structure, which 
clearly is closely connected to the individual eigenvalue 
peaks at the microscopic scale where units are given by 
the average eigenvalue spacing. We thus expect that the 
distribution of the fc-th smallest Dirac operator eigen- 
value corresponds closely to the fc-th peak in the mi- 
croscopic eigenvalue density, a feature that indeed had 
already been noticed in the case of just the smallest 
eigenvalue [23]. In particular, we expect that the fc-th 
eigenvalue distribution as computed in the framework of 
Random Matrix Theory is universal, which indeed turns 
out to be the case (see below). Of course, as one further 
additional check, the sequential summing-up of the in- 
dividual eigenvalue distributions should simply build up 
the full microscopic eigenvalue densities, which are now 
known in closed analytical form. 

The chiral Random Matrix Theory ensembles are de- 
fined by [3-5]: 
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where the integrals are over complex, real, and quaternion real (N + v) x N matrices W for (3 = 2, 1,4, respectively. 
These chiral random matrix ensembles provide exactly equivalent descriptions of the effective field theory partition 
functions in the microscopic finite- volume scaling regime [3,5]. Moreover, all microscopic spectral properties of the 
Dirac operator coincide exactly with the corresponding microscopically rescaled Random Matrix Theory eigenvalues 
[8]. This means that either formulation can be used to derive identical results, and in the case of the individual 
distributions of the smallest Dirac operator eigenvalues it is presently more simple to use the Random Matrix Theory 
formulation. Since the results turn out to be universal, i.e. independent of the detailed form of the Random Matrix 
Theory potential v(WtW) [24-26], it suffices for us to concentrate on Gaussian ensembles with v(x) — x. This choice 
leads to Wigner's semi-circle law p(X) = (2/7r)v2iV — A 2 . 

The partition functions (1) of such chiral Gaussian ensembles corresponding to Nf massive flavors and topological 
charge v can then be written in terms of eigenvalues {a:,} of the positive-definite Wishart matrices W^W, 
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up to an overall irrelevant factor which is independent of the m's. Because everything will be symmetric under 
v — ► —v, we for convenience take v to be non-negative. In the following, we impose the one single technical restriction 
that for j3 = 1 the topological charge v be odd. 

The unnormalized joint probability distributions for N eigenvalues {0 < x\ < ■ ■ ■ < xjy} in the above Random 
Matrix Theory ensembles take the form 

pf(x lt ...,x N ; {m 2 }) = n(«f (V+1) ' 2 - 1 ^ J[{xi + m 2 )) {[ \x t xtf . (3) 
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The joint probability distribution of the k smallest eigenvalues {0 < x\ < ■ ■ ■ < Xk-i < x k } can be written as (the 
order of Xi, . . . , Xk-i can be relaxed): 
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where stands for the normalizing integral 
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s w } ({™ 2 }) = J dxi---J dx N p^\xi,...,x N ;{m 2 }) . (5) 
We now shift Xi — > Xj + X& in the integrand of (4): 
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To get the probability distributions of the Dirac operator eigenvalues we must change the picture back to chiral 
Gaussian ensembles, and take the microscopic limit 
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' T % i — > oo, with the rescaled variables ^ ^(OK/xT kept fixed. (7) 
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In this large- iV limit the difference between partition functions based on N — k and TV eigenvalues becomes insignificant. 
Moreover, one notices that the new terms in the integrand of (6) can be interpreted as arising from new additional 
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fcrmion species, with the partition function now being evaluated in a fixed topological sector of effective charge 
v = 1 + 2/(3. Taking into account the definition (2), we then get, with Z^\{^\) denoting the partition functions in 
the microscopic limit, 
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This is the main result of our paper. It shows that the joint probability distribution of the first fc eigenvalues is given, 
apart from the relatively simple prefactor, by a ratio of finite- volume partition functions. The new partition function 
that enters in the numerator of Eq.(8) has the Nf original fermion masses shifted according to in — > \J ' \i\ + 



contains /3(fc — 1) new fermions of masses — Ci j • • • i y Cfc — Cfc_i (each mass /3-fold degenerate), (3{v + l)/2 — 1 
additional degenerate fermion species of common mass Cfc , and this whole partition function is evaluated in a sector of 
fixed topological charge v = 1 + 2/ f3. While this topological index is fractional for (3 = 4, there is no difficulty with the 
evaluation of the pertinent partition function. (Indeed, it can alternatively be viewed as a partition function in a sector 
of topological charge v = and 3 additional massless fermions.) This expression entirely in terms of the effective field 
theory partition functions strongly suggests that it should be possible to derive these analytical expressions starting 
directly from the effective field theory, perhaps partially quenched as in Rcf. [8]. The proportionality constant in 
Eq.(8) depends on the normalization conventions of the involved partition functions. We fix this numerical factor 
uniquely by the requirement that the total probability of finding any given eigenvalue is normalized to unity. 

Correlation functions of the rescaled eigenvalues {Cki , . . . , £fc p _i , Cfc} are obtained from u)ff\(^i, . . . ,(k', {a 1 }) by inte- 
grating out the remaining eigenvalues in a cell < (i < • • • < Ck ■ In particular, the distribution of the fc-th smallest 
eigenvalue £, p]f\Ci {m})> is given by 
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We note that this gives a very simple representation of the probability distribution of the k-th smallest Dirac operator 
eigenvalue. In particular, already for the most elementary case of just the lowest Dirac operator eigenvalue in the 
(3 = 2 universality class we obtain a very simple expression for an arbitrary topological sector v. In the normalization 
convention of, e.g., [7], 
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a more compact and convenient expression than the one 
provided in Ref. [23]. For example, in the quenched case 
of Nf = it yields the simple relation 

pl 2) (0 = ^e-f 2 /4 d et[7 2+i _ j( C)] , (11) 

where the determinant is over a matrix of size v x v. It is 
worth pointing out that also the corresponding quenched 
distributions for (3 = 1 become exceedingly simple: 

p?\0 = const. C (3-)/2 c -C 2 /8p f W 3 (C) ] , 

(12) 



where the indices i and j run between — vj2 + 1 and 
i>/2 — 1, and the Pfaffian is thus taken over a matrix of 

size [y — 1) x [y — 1) [2]. For example, for v = 1 this 
formula gives 

^(0 = ^/8, (13) 
for v = 3 we obtain 
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while for v = 5 it gives 

(3/ 3 (C) 2 -4/ 2 (C)/4(C) + /i(C)/5(C)) , 

(15) 



all of these being correctly normalized. 

By now, the effective finite- volume partition functions 
are known in analytically closed forms for (3 = 2 [28] 
and (3 — 1 [16]. They are also known for a general even 
number of pairwise degenerate fermions for (3 = 4 [16] 



(this is a technical restriction that we expect to be lifted 
eventually). We shall here for convenience provide the 
precise detailed expressions for all these known cases. In 
exhibiting explicit forms of w[^(£i, . . . , {/i}) by sub- 
stituting these expressions, we choose to set v — 2/ [3 — 1 
in Eq.(2) and introduce additional (3{v+l)/2—l massless 
flavors instead. The number of flavors in this alternative 
picture, Nf + [3(u + l)/2 — 1, will be denoted by n. For 
the f3 = 4 case, we thus consider here only the case with 
an odd number of massless flavors, so that the 'effective' 
number of massless flavors is even, n = 2a. We then get 
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where the matrices and B^' are given by (Q = — Q), jli = x/Cj 2 + M? ) 
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In all cases the normalization constants are fixed uniquely 
by the requirement that the probabilities sum up to unity. 

By construction, the individual distributions 
pf\C'Ap}) sum U P to tnc microscopic spectral density 

pf (CM): 
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To illustrate this, we plot in Fig.l p k ' (() for k — 1, . . . , 4, 

their sum, and ^(f) = \C\[Jq(C) + ^(C))/2 for the 

quenched (Nf = 0) chiral unitary ensemble with v = 0. 
One clearly sees how the microscopical spectral density 
gradually builds up from the individual eigenvalue distri- 
butions. 

We finally turn to the issue of universality. It was 
proven in Refs. [26,27] that the diagonal elements of the 
quaternion kernels K^ 1,4 ^ (C C'i {p}) f° r the orthogonal 
and symplectic ensembles can be constructed from the 
scalar kernel K ^(C {m}) of a unitary ensemble with 
a related weight function. As the scalar kernel in the 
microscopic limit (7) is insensitive to the details of the 
potential v(x) either in the absence [24,25] or in the pres- 
ence of hnite and non-zero masses it [14], so are the corre- 
sponding quaternion kernels. Furthermore, pf{Ci{p}) 
can be expressed in terms of the Fredholm determinant 
ofi^)(C,C';M) [29,30]: 



Here K^f stand for the integral operators with con- 
volving kernels K^\<;,q'; {/x}) over an interval < 
< C The universality of the probability distribu- 
tion pf{C',{p}) i s hence guaranteed. It is nevertheless 
instructive to see how this manifests itself in our explicit 
computation. For a generic Random Matrix Theory po- 
tential the exponential factor e-P( N ~ k ) x k which is pro- 
duced by the shift simply replaced by 
cxp(-(/3/8)(Trp(0)) 2 x k (l + O(l/N))), and thus yields an 

identical factor e _ ^ fc / 8 in the microscopic limit (7) [23]. 
Based on the universality theorems in Refs. [24,26] one 
readily establishes that also the remaining ratio of parti- 
tion functions, and in particular the full expression (8), 
is universal. 
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FIG. 1. Microscopic spectral density ps(C) (thin line), fc-th smallest eigenvalue distribution Pfc(C), k = 1,2,3,4 (solid lines), 



and their sum ^2 k=1 Pk(0 (broken line) for the quenched chiral unitary ensemble (f3 — 2, Nf 



0). 



G 



